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A model of a three Pomeron contribution to high energy elastic pp and pp scat- 
tering is proposed. The data are well described for all momenta (0.01 < \t\ < 
14. GeV 2 ) and energies (8. < «/s < 1800. GeV) (x 2 / d - oi - = 2 -74). The model 
predicts the appearance of two dips in the differential cross-section which will 
be measured at LHC. The parameters of the Pomeron trajectories are: 
a(0) Tl = 1.058, a'(0) Pl = 0.560 (GeV~ 2 ); 
a(0)p 2 = 1.167, a'(0)p 2 = 0.273 (GeV ); 
a(0)p 3 = 1.203, a'(0)p 3 = 0.094 (GeV ). 



1 INTRODUCTION 

Fervently awaited high-energy collisions at LHC will give an access not only 
to yet unexplored small distances but also simultaneously to neither explored 
large distances [0. Future measurements of total and elastic cross-sections 
at LHC P| tightly related to the latter domain naturally stimulate further 
searches for new approaches to diffractive scattering at high energies. 

Recently some models with multi-Pomeron structures were proposed [|3], 
|5[. Some of them 0], [|J use Born amplitudes with two Pomerons as 
single or double poles . Formal violation of the Froissart-Martin bound 
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in some of these models is considered as "practically negligible" though in 
terms of parial-wave amlitudes unitarity violation is flagrant at present-day 
energies. Nonetheless a model of such kind []5| based on the two-Pomeron 
approach shows quite a good agreement with DIS data. 

The eikonal models that are capable of describing the data for nonzero 
transferred momenta are developed in Refs ||, 0. In some cases a "gener- 
alized eikonal representation" is used || together with a dipole (monopole) 
Pomeron contibution, in the others the conventional eikonal is supplemented 
with a "QCD motivated" part consisting of three terms f7j. It is worth notic- 
ing that the two-Pomeron eikonal has been applied to the description of the 
data more than ten years ago (see, e.g., Ref. ||). 

The very multiformity of the models hints that maybe the most general 
way to describe high-energy diffraction is just to admit an arbitrary num- 
ber of Pomerons (i.e. all vacuum Regge-poles contributing non-negligibly at 
reasonably high energies. Roughly, they should have intercepts not lower 
than 1). On the one hand this seems not very economical. But on the other 
hand we could argue that no basic principle forbids more than one single 
Pomeron. We could also add that in the perturbative framework the account 
of the renormalization group leads presumably to converting of the fixed 
branch point (in the J-plane) into an infinite series of simple poles accumu- 
lating down to 1 from some maximal value ||. Unfortunately perturbative 
searches in this field are far from being satisfactory from many viewpoints. 

In this paper we would like to make a first step in realization of the above 
formulated hypothesis about many-Pomeron structure of the eikonal. As 
it seems not possible to describe the data in the framework of the eikonal 
approach with presence of one single pole Pomeron contribution ||10|| , and 
the two-Pomeron option does not improve quality of description drastically 
(more details are given in the text) it is fairly natural to try the next, three- 
Pomeron, option for the eikonal. We will see below that this choice appears 
rather lucky. 

2 THE MODEL 

Let us brifely outline the basic properties of our model. Unitarity condition: 

3mT(s,6) = \T(s,b)\ 2 + r](s,b) , 
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where T(s, b) is the scattering amplitude in the impact representation, b is 
the impact parameter, r](s, b) is the contribution of inelastic channels, implies 
the following eikonal form for the scattering amplitude T(s, b) 

2iS(s,b) _ 1 

T(s,b) = , (1) 

— # 

where S(s, b) is the eikonal function. The unitarity condition in terms of the 
eikonal looks as follows 

Sm 5{s, b) > 0, s > s ine i . (2) 

The eikonal function is assumed to have simple poles in the complex J- 
plane and the corresponding Regge trajectories are normally being used in 
the linear approximation 

a(t) = a(0) + a'(0)t . (3) 

Accordingly we get the following (modulo the signature factor) contribu- 
tion to the eikonal function in t-space (here t is the momentum transfer) 

where 

p 2 = 4a'(0)ln- + r 2 (5) 

So 

is referred to as the "Reggeon radius" . 

In order to relate t- and 6-spaces one proceeds via Fourier-Bessel trans- 
forms 

/■oo 

f(t)=4ns db 2 J (bV^i)f(b) , 



1 
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(6) 
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Making use of Eq. ([|) we obtain the following 6-representation of the 
eikonal function M) 



For the cross-sections we use the following normalizations: 



otot = -$smT(s,t = 0), 
s 



POO 

(^elastic = 4vr / db 2 \T(s,b)\ 2 , 
Jo 

(8) 



da \T(s,t)\' 



dt lQns 2 

9tfeT(s,t = 0) 



P 



SsmT(s,t = 0) 



In the present model we assume the following representation for the 
eikonal function: 

5%{s, b) = 8+ (s, b) + 5+ (s, b) + 5+(s, b) T 5®(s, b) + 8+(s, b) =f <£(s, b), (9) 

here 5p (s, 6) are Pomeron contributions. '+' denotes C even trajecto- 
ries (the Pomeron trajectories have the following quantum numbers + J ++ ), 
'— ' denotes C odd trajectories, Sq(s, b) is the Odderon contribution (the Odd- 
eron is the C odd partner of the Pomeron with quantum numbers 0~J ); 
Sj', 6~(s,b) are the contributions of secondary Reggeons, f (C = +1) and uj 
(C = -l). 

The form (|]) is not compatible with analyticity and crossing symmetry, 
which are easily restored by substitution s — > se~ l7T ^ 2 . We introduce a new 
dimensionless variable 

s = -e^i , (10) 
So 

and obtain each C+ and C— contribution with its appropriate signature 
factor and the form: 



5 + (s,b) = i 



£ja(0)-l£_ 



s Airp 2 
p 2 = 4a'(0)lnS + r 2 , 
(C = +l); 



(11) 
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r(s,b) 
P 



_b 2 
2 ' 



s 4vrp 
2 - 4a'(0)lns" + r 2 



(12) 



(C = -l). 

The parameters of secondary Reggeon trajectories are fixed according to 



(13) 



the parameters obtained from a fit of the meson spectrum |TT 

a f (t) = 0.69 + 0.84t , 

a u (t) = 0.47 + 0.93t . 

All the trajectories are taken in linear aproximation: 

ai (t) = ai (0) + a[(0)t, (i = P 1 ,P 2 ,P 3 ,0). (14) 

Let us remark that all good fits require ap(0) — 1 = Ap > which means 
that Born amplitude will eventually exceed the Froissart-Martin ]T2| unitar- 
ity bound. This violation of unitarity is removed by all kinds of "eikonal- 
ization". Nevertheless, one must take into account the following unitarity 



constrains 13 



a P (0) > a o (0) and a' ¥ (0) > a' o {0) , (15) 

there P in this case is the leading Pomeron trajectory (the one with the 
highest intercept A F ). 



3 RESULTS 

We fitted the adjustable parameters over a set of 961 pp and pp data of 
both forward observables (total cross-sections a ta t, and p - ratios of real to 
imaginary part of the amplitude) in the range 8. < \/s < 1800. GeV and 
angular distributions (^) in the ranges 23. < sfs < 1800. GeV, 0.01 < \t\ < 
14. GeV 2 . 

Having used 20 adjustable parameters we achieved xVd.o.f. = 2.74. The 
parameters are presented in Table [I] (all the errors are obtained according to 
MINUIT output). 
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Pomerorii 




f-Reggeon 




0.0578 ±0.0020 


A f 

J 


-0.31 (FIXED) 


QPi 


53.007 ±0.795 


c f 

J 


191.69 ±2.12 


aL 

iri 

< TIE 


0.5596 ±0.0078 (GeV' 2 ) 
6.3096 ±0.2522 (GeV' 2 ) 


i 
r 2 f 


0.84 (GeV- 2 ) (FIXED) 
31.593 ± 1.099 (GeV~ 2 ) 




Pomerori2 




cu-Reggeon 




0.1669 ±0.0012 


A w 


-0.53 (FIXED) 




9.6762 ±0.1600 




-174.18 ± 2.72 


a™ 


0.2733 ±0.0056 (GeV- 2 ) 


a' 

Ul 


0.93 (GeV- 2 ) (FIXED) 
7.467 ± 1.083 (GeV' 2 ) 


^2 

/ TIE 

Jr2 


3.1097 ±0.1817 (GeV- 2 ) 


r 2 , 




Pomeron 3 






A P3 


0.2032 ±0.0041 


so 


1.0 (GeV 2 ) (FIXED) 


Cp 3 


1.6654 ±0.0669 






«p 3 


0.0937 ±0.0029 (GeV- 2 ) 








2.4771 ± 0.0964 (GeV- 2 ) 








Odderon 






A 


0.19200 ±0.0025 






Co 


0.0166 ±0.0022 






«o 


0.048 ± 0.0027 (GeV- 2 ) 






r 2 


0.1398 ±0.0570 (GeV~ 2 ) 







Table 1: Parameters obtained by fitting to the data. 

In order to estimate the quality of the description, we have calculated 
partial \ 2 over an sets of data used in the fit. This x 2 is calculated using the 
following formula: 

2 _ \^ (Vtheory(ri) ~ 0~exp(ri)) 2 . , 

X ~tl (A(a exp (n))) 2 ' 1 } 

where ntot is the number of data in the set, a exp is the experimental value of 
the quantity that is described, o~theor y is our prediction for this quantity, and 
A(cr ea; p(n)) is the experimental uncertainty. 
The partial \ 2 may be found in Table 0. 
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Set of data 


Number of points, ntot 


X 2 / ntot 


1 

2 
3 
4 
5 
6 
7 
8 
9 

10 
11 
12 
13 
14 
15 
16 


pp 

u total 

PP 
u total 
pPP 

pPP 

v/i = 31. (GeV) 
d £\ v/i = 53. (GeV) 
§ pp , yCs = 62. (GeV) 
s| pp , = 546. (GeV) 
= 630. (GeV) 

^ = 1800. (GeV) 
§ PP , = 23.5 (Gel/) 
§ pp , ^ = 27.43 (GeV) 
f pp , v^ = 30.7 (GeV) 
*L m t = 44.64 (GeV) 
= 52.8 (GeV) 

^ = 62. (Gey) 


42 

50 

11 

36 

22 

52 

23 

78 

19 

51 

105 

39 

92 

97 

93 

151 


2.3524 
0.6309 
0.6942 
1.9075 
3.3688 
8.5457 
1.8524 
3.8425 
9.9273 
1.3741 
2.2491 
1.8929 
4.4559 
1.5748 
2.0956 
2.4272 




Number of parameters 


Total number of points 


xVd.o.f. 




20 


961 


2.7441 



Table 2: Partial \ 2 ■ 



Some of these x 2 s are high (for instance those for differential cross sections 
at \fs = 53, 630 GeV). It reflects the fact that we did not make use of 
systematical errors for these sets of data which can be as high as 30%. 

The results are shown in fig. |, |, § (§ §, § 0, g 

We do not include elastic cross-section data sets into the fit and predic- 
tions of the model for elastic cross-sections can be seen in fig. 

It is instructive to compare these results with a two-Pomeron option. We 
give corresponding results in fig. |9], [H], [TT], [l^, [13|. In this case x 2 /d.o.f. = 
10.87 what can be considered as an argument in favour of our hypothesis 
of many Pomerons. The parameters obtained in the two-Pomeron option 
are presented in Table |^ (all the errors are obtained according to MINUIT 
output). 
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Pomeroni 




f-Reggeon 


QPi 


0.0859 ±0.0021 
53.18 ±0.86 


J 

c f 


-0.31 (FIXED) 
188.51 ± 12.13 


aL 

Fx 


0.360 ± 0.009 (GeV~ 2 ) 
9.595 ± 0.6289 (GeV~ 2 ) 


a' f 

j 

rj 


0.84 (GeV' 2 ) (FIXED) 
41.424 ± 7.971 (GeV~ 2 ) 




Pomeron2 




w-Reggeon 


Ap 2 


0.14437 ±0.0051 


A w 


-0.53 (FIXED) 


cp 2 


6.87 ±0.36 




-171.36 ±8.23 


«p 2 

y < 2 

'P 2 


0.082 ± 0.004 (GeV 2 ) 
4.765 ± 0.2533 (GeV- 2 ) 


< 


0.93 (GeV- 2 ) (FIXED) 
2.621 ±6.362 (GeV' 2 ) 




Odderon 






A 

Co 


-0.2707 ±0.1178 
1.8134 ± 1.4837 


so 


1.0 (GeV 2 ) (FIXED) 


' O 


0.029 ± 0.023 (GeV' 2 ) 
1.159 ±0.591 (GeV- 2 ) 







Table 3: Parameters obtained by fitting to the data with two Pomeron con- 
tributions. 
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Figure 1: Total cross sections of pp scattering (hollow circles) and pp scatter- 
ing (full circles) and curves corresponding to their description in the present 
model. 



9 















} 












Jfl 

n 








1 / 




If 


f 

ill 
If 


i. 


P-P-- 




• / 


LHC 


















It 












\n 


r 1 1 i jrfr 












l^------------- 










P-P-- - - 








1 













1 10 1 2 1 3 10* 



Vs (GeV) 

Figure 2: Elastic cross sections of pp scattering (hollow circles) and pp scatter- 
ing (full circles) and curves corresponding to their description in the present 
model. These sets of data are not included in the fit. 
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Figure 3: Ratios of the real to the imaginary part of the forward pp scatter- 
ing amplitude (hollow circles) and pp scattering amplitude (full circles) and 
curves corresponding to their description in the present model. 
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-t (GeV 2 ) 

Figure 4: Differential cross-sections for pp scattering and curves correspond- 
ing to their description in the present model. A 10~ 2 factor between each 
successive set of data is omitted. 
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Figure 5: Differential cross-sections for pp scattering in the region of small 
momenta 0.01 < \t\ < 1.2 GeV 2 and curves corresponding to their description 
in present the model. A 10~ 2 factor between each successive set of data is 
omitted. 
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Figure 6: Differential cross-sections for pp scattering and curves correspond- 
ing to their description in the present model. A 10~ 2 factor between each 
successive set of data is omitted. 
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Figure 7: Differential cross-sections for pp scattering in the region of small 
momenta 0.01 < \t\ < 1.2 GeV 2 and curves corresponding to their description 
in the present model. A 10~ 2 factor between each successive set of data is 
omitted. 
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Figure 8: Predictions of the model for the differential cross-section of pp 
scattering which will be mesured at LHC with y/s = 14. TeV and at RH1C 
y/s = 500. GeV . The data corresponding to the energy y/s = 62. GeV is 
multiplied by 10" 8 , RHIC by 10~ 12 , and that of LHC by 10~ 16 . 
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Vs (GeV) 

Figure 9: Total cross sections of pp scattering (hollow circles) and pp scat- 
tering (full circles) and curves corresponding to their description in the two- 
Pomeron model. 
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Figure 10: Ratios of the real to the imaginary part of the forward pp 
scattering amplitude (hollow circles) and pp scattering amplitude (full circles) 
and curves corresponding to their description in the two-Pomeron model. 
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Figure 11: Differential cross-sections for pp scattering and curves correspond- 
ing to their description in the two-Pomeron model. A 1CT 2 factor between 
each successive set of data is omitted. 
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Figure 12: Differential cross-sections for pp scattering and curves correspond- 
ing to their description in the two-Pomeron model. A 10~ 2 factor between 
each successive set of data is omitted. 
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-t (GeV 2 ) 

Figure 13: Predictions of the two-Pomeron model for the differential cross- 
section of pp scattering which will be mesured at LHC with y/s = 14. TeV 
and at RHIC y/s = 500. GeV. The data corresponding to the energy y/s = 
62. GeV is multiplied by 10" 8 , RHIC by 10~ 12 , and that of LHC by 10~ 16 . 



19 



4 CONCLUSION AND DISSCUSSION 



Above we have developed a model which is based on a general argument of 
multiplicity of the Pomeron Regge poles in the eikonal. The present model 
shows a very good descriprion of the available data for all momenta (0.01 < 
\t\ < 14. GeV 2 ) and energies (8. < < 1800. GeV) so that x 2 /d.o.f. = 
2.74. 

The model predicts the appearance of two dips in the differential cross- 
section which will be measured at LHC, fig. |8|, and this prediction is stable 
in the sence that all two models with two and three Pomeron contributions 
predict the same behaviour of the differential cross-section with two dips. 
These dips are to apear in the region t\ ~ —0.5 GeV 2 and t 2 — —2.5 GeV 2 
which is in agreement with other predictions (model ||). 

We predict the following values of the total cross-section, elastic cross- 
section, and the ratio of real to imaginary part of the amplitude for the 
LHC: 



= 14. TeV 

aZ = 106.73 (mb) ™ ™ b , (17) 



a vp = 9Q 19 (mb) + 3 - 58mb 

elastic ^.LV {"M) _ 2 .83 mb ' 



«* = 0.1378 ±S:88g 



Predictions for RHIC are: 



= 500. GeV , 
og = 59.05 (mb)tlitZl (18) 
12-29 (mb)t° ?^l 



ffP = 0.1327 ±°™ 52 



1.0071 
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The parameters of the Pomeron trajectories are: 



a(0) Pl = 1.058, a'(0) Pl = 0.560 (Gel/" 2 ); 

a(0) P2 = 1.167, a'(0)p 2 = 0.273 {GeV~ 2 )] (19) 

a(0)p 3 = 1.203, a'(0) P3 = 0.094 (GeV~ 2 ). 

Their coupling constants fulfil the following inequality: 

c Pl > cp 2 > cp 3 (20) 

It may mean that there exists a series of Pomeron contributions each 
term having the form (JXTJ), or these three contributions effectively emulate 
one nonlinear Pomeron trajectory. 

The intercepts and slopes fulfil the following inequalities: 

A Pl < A P2 < Ap 3 (21) 

a Pl (o) > a P2 (o) > 4 3 (o), 

i.e. the higher is the intercept the lower is the slope. We observe that the 
product of the intercept and the slope is approximately the same for all the 
Pomerons, A ■ a'(0) ~ 0.040 ± 0.0009 (GeV~ 2 ). This is seen in fig. |TJ. This 
constant seems suprisingly universal if compared with the products of other 
Reggeon parameters used in this model (fig. [r5| ). At present we have no clear 
understanding of this universality. 

We can only remind that high energy asymptotic behaviour of total and 
elastic cross-sections in Regge-eikonal approach have the following form: 

a tot (s) ^87ra P (0)A P ln 2 (s/s ) , (22) 

CTelastic(s) -> 47TCtp(0) A P \b. 2 (s/s ) , 

and the constant a P (0)Ap (GeV~ 2 ) (there P stands for the rightmost singu- 
larity of the eikonal function in J-plane) defines a universal (independent on 
colliding beams) asymptotic behaviour. 
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Figure 14: Products of intercepts 
to slopes for the Pomerons (hol- 
low squares) and the Pomeron in 
a generalized eikonalization model 
Pdga/p [fj, Donnachie&Landshoff 
Supercritical Pomeron (hol- 
low triangles). Solid line corre- 
sponds to the mean value of the 
product (for the three Pomerons) 
and dashed lines correspond to its 
error corridor. 



Figure 15: Products of intercepts 
to slopes for the Pomerons, the 
Odderon, and Reggeons used in 
the present model. 



It is interesting to enlist the following characteristic properties of the 
Pomerons used in this paper. 

The first of the Pomerons ('Pomeroni) possesses the properties that we 
expect from the string picture ]n| of Reggeons, i.e. a'(0)p = |a'(0)/- = 
0.42 (GeV~ 2 ) and indeed a'(0) Pl = 0.559 ± 0.078 (Gel/- 2 ). 

The second Pomeron ('Pomeron^) is close to what is called "supercritical 
Pomeron" with the slope a'(0) Pa = 0.273±0.005 (GeV~ 2 ) close to its "world" 
value a'(0) P ~ 0.25 (Gel/- 2 ). 

The third Pomeron ('Pomerons) is reminiscent of what is known as a 
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"hard" ( or perturbative QCD )Pomeron. Its parameters (a(0)p 3 = 1.203, 
a'(0)p 3 = 0.094 (GeV~ 2 )) are close to the calculated parameters of the per- 
turbative Pomeron, which arise from the summation of reggeized gluon lad- 
ders and BFKL equation @: a(0)^ FKL ~ 1.2, a'(0)^ FKL ~ 0. (GeV~ 2 ). 
The fact of arising of a "hard" Pomeron in a presumably "soft" framework 
can seem quite unexpected. However we are not particularly inclined to iden- 
tify straightforwardly "our hard Pomeron" with that which is a subject of 
perturbative QCD studies. 

The Odderon has the following parameters: a(0)® = 1.192, a'(0)o = 
0.048 (GeV~ 2 ) in agreement with unitarity constraints Eq. (0). The Odd- 
eron intercept is positive and close to that of the Pomeron?,. The slope is 
almost zero. The coupling is so small that only high-t data may be sensible 
to the Odderon contribution. 

Assuming that one can neglect the non-linearities of Regge trajectories 
and making use of a simple parametrization 

a(m 2 ) = a(0) + a'(0) ■ m 2 , (23) 

we can try to estimate the corresponding spectroscopic content of our model. 

Then 3ftea(m 2 ) = J, where J is an integer number corresponding to the 
spin of a particle which we should find lying on the trajectory. 

The trajectories are depicted in fig. 0. The C+ Reggeon trajectory is 
in fact a combination of two families of mesons / and a 2 . The C— Reggeon 
trajectory is a combination of two families of mesons uj and p. As is seen, the 
secondary Reggeon trajectories fairly well describe the spectrum of mesons. 

Among the mesons with appropriate quantum numbers there exesit two 
that fit the Pomeron trajectory (0 + J ++ ): /2(1810) + 2 ++ with mass m = 
1815 ±12 MeV and X(1900) 0+2++ with mass m = 1926 ± 12 MeV . One 
of them is supposed to be on Pomeron 2 trajectory. 
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Figure 16: Regge trajectories of secondary Reggeons, 
three Pomerons and the Odderon. 
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